FINITE-TIME BLOWUP FOR THE INVISCID PRIMITIVE EQUATIONS OF 
OCEANIC AND ATMOSPHERIC DYNAMICS 

CHONGSHENG CAO, SLIM IBRAHIM, KENJI NAKANISHI, AND EDRISS S. TITI 

Abstract. In an earlier work we have shown the global (for all initial data and all time) well-posedness of 
strong solutions to the three-dimensional viscous primitive equations of large scale oceanic and atmospheric 
dynamics. In this paper we show that for certain class of initial data the corresponding smooth solutions of 
the inviscid (non-viscous) primitive equations blow up in finite time. Specifically, we consider the three- 
dimensional inviscid primitive equations in a three-dimensional infinite horizontal channel, subject to 
periodic boundary conditions in the horizontal directions, and with no-normal flow boundary conditions 
on the solid, top and bottom, boundaries. For certain class of initial data we reduce this system into 
the two-dimensional system of primitive equations in an infinite horizontal strip with the same type of 
boundary conditions; and then show that for specific sub-class of initial data the corresponding smooth 
solutions of the reduced inviscid two-dimensional system develop singularities in finite time. 
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1. Introduction 

The three-dimensional primitive equations for large scale oceanic and atmospheric dynamics are given 



by the system of partial differential equations: 

u t + uu x + VUy + wu z + p x - Rv = v H A H u + v 3 u zz , (1) 

v t + uv x + Wy + wv z + py + Ru = v H A H v + v 3 v zz , (2) 

Pz + T = 0, (3) 

T t + uT x + vT y + wT z = Q + k h A h T + k 3 T zz , (4) 

u x + v y + w z = 0, (5) 



which are supplemented with the initial value (uo,Vo,Tq), and satisfy the relevant geophysical boundary 
conditions (see, e.g., [16, 19, 20]. Here Ah = d xx + d yy denotes the horizontal Laplacian operator. 
The global well-posedness (for all time and for all initial data) of strong solutions to the above three- 
dimensional system, subject to the relevant geophysical boundary conditions, has been proven first in [6] 
with full viscosity (i.e. vh > and v 3 > 0) and full diffusion (i.e. kjj > and n 3 > 0), (see also, [14] 
for the case of Dirichlet boundary conditions). This result has been improved recently in [7] for the case 
of full viscosity (i.e. vh > and v 3 > 0) and only partial anisotropic vertical diffusion (i.e. kh = and 
k 3 > 0) which stands for the vertical eddy heat diffusivity turbulence mixing coefficient (see, e.g., [12], [13], 
for the geophysical justification). In the above results, the Coriolis forcing term, with rotation parameter 
R, did not play any role in proving the global regularity. This is contrary to the cases of the three- 
dimensional fast rotating Euler, Navier-Stokes and Boussinesq equations by [2, 3, 4, 5] where the authors 
take full advantage of the absence of resonances between the fast rotation and the nonlinear advection. 
This absence of resonances at the limit of fast rotation leads to strong dispersion and averaging mechanism 
(at the limit of fast rotation for large values R depending on the size of the initial data) that weakens the 
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nonlinear effects and hence allows for establishing the global regularity result in the viscous Navier-Stokes 
case, and prolongs the life-space of the solution in the Euler case (see also [8, 11] and references therein; 
in addition, see [1] for simple examples demonstrating the above mechanism). 

In geophysical situations it is observed that the viscosity coefficients are very small, and that in fact 
f3 C « 1. Motivated by this observation and the above discussion it will become interesting to know of 
whether the inviscid (non-viscous) primitive equations are globally regular or that they develop singularity 
(blow up) in finite time. Since the rotation term did not play any role in establishing the global regularity 
in the viscous cases (cf. [6, 7]) one might as well ask the above question without the rotation term (i.e., 
consider the case with R = 0). Therefore, we will consider in this paper the inviscid primitive equations 
without the Coriolis rotation term, and we will show that for certain class of smooth initial data their 
corresponding smooth solutions will develop a singularity (blowup) in finite time. 

In order to establish our result we will assume that we are given a smooth solution to the inviscid 
primitive equation. In section 2 we will derive a reduced equation that this smooth solution will satisfy. 
In section 3 we will follow [9] (see also [17]) to show that for certain class of initial data the corresponding 
solutions to this reduced equation blow up in finite time. In section 4 we provide a family of initial data 
whose corresponding smooth solutions to the inviscid primitive equations blow up in finite time. It is 
worth mentioning that similar approach has been introduced in [10] to show the blowup for the Prandtl 
equation of the boundary layer in the upper half-space. 



2. Derivation of a reduced equation 

In this section we will assume that we are given a conveniently smooth solution to the inviscid primitive 
equations. We will derive a reduced equation that this solution must satisfy by requiring that the solution 
fulfills certain symmetry conditions. Eventually, in section 3, we will show that for certain class of initial 
data the corresponding solutions of this reduced equation develop singularity in finite time. 

First, let us consider the inviscid primitive equations without the Coriolis force: 

u t + uu x + vu y + wu z + p x = 0, (6) 

v t + uv x + vv y + wv z + p y = 0, (7) 

Pz +T = 0, (8) 

T t + uT x + vT y + wT z = k h A h T + k 3 T zz , (9) 

u x + v y + w z = 0, (10) 

in the horizontal channel = {(x, y, z) : < z < H, (x, y) 6 M 2 }; subject to the boundary conditions: 
no-normal flow and no heat flux in the vertical direction at the physical solid (top z — H and bottom 
z = 0) boundaries, and periodic boundary conditions, say with period L, in horizontal directions. Observe 
that when the initial temperature T = 0, then it is easy to see that any smooth solution to system (6)- 
(10), subject the above boundary conditions, must satisfy T(x,y,z,t) = 0. Consequently, for any smooth 
solution with T = the velocity field satisfies the following system: 

u t + uu x + vu y + wu z + p x = 0, (11) 

v t + uv x + vv y + wv z +p y = 0, (12) 

Pz = 0, (13) 

u x + v y + w z — . (14) 

Since our goal is to establish the blowup for certain class of smooth solutions and initial data, we will further 
simplify matters and restrict ourselves to smooth solution to the above system that are independent of the 
y— variable and that the y— component of initial velocity vector filed vq = 0. Once again it is clear that 
any smooth unique solution to the above system with these properties will satisfy v(x, z, t) = 0, and that 
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u, w arc functions of the spatial variables (x, z) and of the time t only. Consequently, we restrict ourselves 
further and only consider solutions that satisfy the following two-dimensional sub-system: 

u t +uu x + wu z +p x =0, (15) 
Pz = 0, (16) 
u x + w z =Q, (17) 

in the two-dimensional channel 

M = {(x,z) : < z < H ,x G M}, 

subject to the no-normal flow on the vertical direction, and periodic boundary conditions in horizontal 
x— direction, namely, 

w\ z =h = w \ z =o = °> ( 18 ) 

and 

u(x + L, z, t) = u{x, z, t); p{x + L, z, t) — p(x, z, t); w(x + L, z, t) = w(x, z, t). (19) 

Observe that the space of periodic functions with respect to x with the following symmetry 

u(x, z, t) = — u{— x, z, t);p(x, z, t) = p(—x, z, t); w{x, z, t) — w{—x 1 z, t), (20) 

is invariant under the solution operator of system (15)-(17) subject to the boundary condition (18) and 
(19). Therefore, from now on we will restrict our discussion to the solutions of system (15)— (19) that 
satisfy the symmetry condition (20). 
Let us denote by 

4fat) = jjJ <p{x,z,t)dz. (21) 
From equation (17) and the boundary condition (18), we have 

= w(x, H, t) — w(x, 0, t) = / w z (x, z, t) dz = — f u x (x,z,t)dz. (22) 



Thus, 

u x (x,t)=0. (23) 
By differentiating equation (15), with respect to x, we obtain 

du x 

— + {uu x ) x + w x u z + wu zx +p xx = 0. (24) 

Equation (16) implies that p(x, z,t) is only function of (x,t), i.e., p(x, z,t) — p(x,t). By averaging (24) 
with respect to the z variable over [0, if] and using (23), we obtain 



[uu x ) x + w x u z + wu xz +p xx = 0. (25) 

Notice that (18) implies 

w x (x,H,t)=w x (x,0,t)=0. (26) 
Thanks to (17) and (26), integrating by parts gives 



w x u z + um xz = —w xz u - w z u x 



= u xx u + (u x ) 2 = {uu x ) x . 
As a result, (25) can be reduced to 



2(uu x ) x +p xx =0. (27) 
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Thus, 

2uu^ + p x = C(t), (28) 

for some function C(t). Since p x and uu^ are odd functions, with respect to the variable x (see condition 
(20)), then 

uu^(0,t) =p x (0,t) = 0. 

Therefore, 

C(t) = 2uu^(0,t)+p x (Q,t) = 0, (29) 

and consequently, 

p x = -2uu x ~. (30) 
Substituting (30) into system (15)-(17) we obtain the closed system 

du 

— + uu x + wu z - 2uu x = 0, (31) 
at 

u x + w z =0, (32) 

subject to the boundary conditions (18)-(20). In particular, by differentiating with respect to x we have 

flu 

+ (uu x ) x + w x u z + wu xz - 2(uu x ) x = 0. (33) 

Let us consider the restriction of the evolution of equation (33) on the line x = 0. Since u(x,z,t) is an 
odd-function, and w(x,z,t) is an even-function, with respect to the variable x we have: 

u(0,z,t) = 0;w x (0,z,t) = 0. 

This together with (33) imply 



u tx (0, z, t) + (u x (0, z, t)f + w(0, z, t)u xz (0, z, t) ~ 2(u x (0, z, t)f = 0. (34) 
Recalling that w z (0, z, t) = — u x (0, z, t), and denoting by W(z, t) — u>(0, z, t), we obtain 

W tz -(W z ) 2 + WW zz + ^ (W z ) 2 dz = 0, (35) 
with the boundary conditions: 

W(0,t) = W(H,t) = 0. (36) 
This is a closed equation that we will investigate in section 3. 

3. Self-similar blowup of the reduced equation 

In this section wc will investigate the blowup of the nonlinear integro-differential boundary value problem 
(35)-(36). The analytical and numerical investigation of this type of equations have been the subject matter 
of the papers [17, 18]; see also [15] for additional new results. The presentation below follows [9] (see also 
[17] section 4, and references therein), where the same problem, arising in a different fluid dynamical 
context, has been investigated. For the sake of completeness we provide below the full details. 

Noting that equation (35) is invariant for the scaling 

W(z,t) i y \W(z,\t), 
therefore, we look for a self-similar solution in the form 

W{z,t) = with^O) = <p{H) = 0. 
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This solution starts from the initial profile W(z, 0) = <p(z), and blows up at time t = 1 at every z G (0, H), 
for which <p(z) ^ 0. Therefore, the idea is to construct a nontrivial initial profile p(z). Using the above 
ansatz and substituting in (35) we obtain the following reduced nonlinear nonlocal boundary value problem: 

iff - W f + W + I j H {^'{z) fdz = 0, ¥>(<)) = p(H) = , (37) 

that we need to show that it has a nontrivial solution ip. 

Let m > be a given free parameter, we consider instead of (37) the nonlinear boundary value problem 

iff - {ip') 2 + ipip" + to 2 = 0, ip(0) = ip(H) = 0. (38) 

We observe that in case the boundary value problem (38) has a solution, which we denote by tp m , then 
ip m is nontrivial because m > 0. Moreover, integrating equation (38) over < z < H yields 

m 2 = \ j" y m {z)fdz. 

Consequently, p m is also a nontrivial solution to the original nonlocal boundary value problem (37). 
Therefore, we will focus now on showing that problem (38) has a nontrivial solution for every m > given. 

First we look for a positive solution ip to the boundary value problem (38). Denoting by ip := <p', 
then any local smooth solution to the second order ordinary differential equation given in (38) can be 
parameterized by an integral curve (ip(z),tp(z)), in the phase portrait space (ip,tp). Suppose that we 
can parameterize, say locally, the integral curve in terms of ip instead of z, i.e., that we can invert the 
relationship between ip and z and express z as function ip. Then equation (38) becomes 

1 dp ip ip 

ip dip ip 2 — ip — m 2 (ip — ip+)(ip — ip~) ' 

where ip± are the singular points 

ip+(m) := y/m 2 + 1/4 + 1/2 > 0, ip-(m) := -yjm 2 + 1/4 + 1/2 < 0. 
Integrating the above equation we obtain general integral curves in the phase portrait for (tp, ip) 

cp = C\ip- ip+\ \ip-ip-\ , ( 39 ) 

with an integration constant C . Since we are interested in positive solutions <p we take C > 0. In order 
for this curve to be a solution to the boundary value problem (38), and that it can be parameterized in 
terms of z G [0, H], one has to require that for each C > the above curve yields a solution satisfying 

(p,iP)(z = 0) = (0,iP+), (p,iP)(z = H) = (0,iP_). (40) 
Injecting (39) into (38) yields 

dib —1 ~*- ^+ 

— = — \ip -iP + \^+~^- \ip - ip_\^+-^- , (41) 
dz C 

which shows that ip is a decreasing function and hence the inversion between the variables z and ip is valid 
over the interval ip_ < ip < ip + . In particular, for every ip G [ip-,ip + ] and by virtue of (41) and (40) we 
have 

m= 77#= / -y— 4 2# = -C/ \ijj-iP + \^-\iP-ip.\^-dip. (42) 

Jtb+ a W Jtb+ V -rV> —m Jib, 
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Moreover, from the other boundary condition in (40) the length of interval is determined by 

f<P+ f- —4>+ 

= C \ip - ip+\ ^+"^- \ip-ip-\ '* + ^- dip 
c _ B ( V+ -<P- 



(43) 



ip + -ip_ \ip + -ip_ ip+-ip- 

C Jl. 1 1 



2y/m 2 + 1/4 \2 A^Jm 2 + 1/4' 2 4 v / m 2 + 1/4 J ' 

where B(-, •) is the Beta function. 

Hence with the given H > 0, for each to > there is a unique C > 0, depending on _ff and m, such 
that (43) holds. 

From all the above we concluded that there is a solution ip of (38) satisfying 

,-H 



2 f 

ip(0)=0 = <p(H), ip>0on0<z<H, and — / (ip'(z)) 2 dz = to 2 . 

H Jo 



This in turn also shows that ip is a nontrivial solution to the nonlinear nonlocal boundary value problem 
(37). 

From (41), (40) and the above discussion one also concludes that 

<Pzz(0)=0 = <p xz (H). (44) 

Furthermore, the symmetry of equation (38) implies that —<fi(H — z) is a negative solution to (38). 

In conclusion, we have obtained a one-parameter family of blowup solutions of (35)-(36), which blow 
up at every z £ (0,H), as t — > 1~. 

Remark 1. If we choose the constant C in (43) such that z(ip-) — 4f , then condition (44) holds at z = ^ 
instead of z = H. This in turn will allow us to glue the positive solution defined on the interval [0, 
and which is also satisfying (40) z = ^ instead of z = H, with its negative counterpart at z = If. Hence, 
this idea will allow us construct sign-changing blowup solutions. In fact, this is the type of solutions that 
are constructed in [9] and [17]. Of course one can repeat the above and get profiles that change signs as 
many times as one wants. 

4. Blowup of the smooth solutions to the primitive equations 

In this section we will demonstrate our main result and show that for certain class of initial data the 
corresponding smooth solutions to the primitive equations (15)-(17), subject to the boundary conditions 
(18)-(20), blow up in finite time. First, let us prove the following proposition concerning the uniqueness 
of solutions to equation (35)-(36). 

Proposition 1. Let Wi,W 2 be two solutions of (35)-(36) which are in L 2 ((0, T); H 2 ), with Wi(z,0) = 
W2(z,0). Then W\(z,t) — W2(z,t), for all t e [0,T). In particular, W(z,t) — ^-j is the only solution of 
(35)-(36) in the space L 2 ((0,T); H 2 ), for allT £ [0, 1), with initial data (fi(z), where ip(z) is any nontrivial 
solution of the boundary value problem (37) that was established in section 3. 

Proof. Let V = W± — W 2l and W = \{Wi + W 2 ). Then (35)-(36) imply: 

V tz - 2W Z V Z + VW ZZ + WV ZZ + ^ J" W Z V Z dz - 0, (45) 
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with the boundary conditions: 

V(0, t) = V(H, t) = 0. (46) 

Multiplying (45), integrating with respect to z over [0, if], integrating by parts and using the boundary 
conditions (36) and (46) gives: 

\jt I" y2{Zi t)dz = \j" Wz{ y^ 2dz - j" w ** vv * dz ■ 

Thus, by Holder inequality we obtain 

2(0, H) ^ 5|| W 2 ||L oo (0,-ff) II (0 IIl 2 (0,-H") + ^\W\\H 2 {0,H)\\y{t)\\L^{0,H)\\yz{t)\\L 2 {0,H) ■ 

By virtue of the one-dimensional Sobolev imbedding theorem, and thanks to the boundary conditions (46) 
one can apply the Poincare inequality, to obtain: 

J t \\VM\\l H 0M) < 4W\\H^,H)\\V z (t)\\lH0M) ■ 

Applying Gronwall's inequality we then conclude that || V z (-, i)) || jl 2 (o ,h) = 0, for all t e ]0,T). Again, 
thanks to (46) we apply Poincare inequality to conclude the uniqueness part of the statement of the 
proposition. The second part of the statement is an immediate corollary of the first part. 

□ 

Theorem 2. Let <p(z) be any nontrivial solution of the boundary value problem (37), and let k be any 
positive integer. Suppose that (u(x, z, t),w(x, z, t)) is a smooth solution of (15)- (11), subject to the boundary 
conditions (18)-(20), with initial value 

Uo(x, z) — -^—^ sm (-^-" a; )¥''( z ) an d u>o(x, z) — cos(^^- x)ip(z) . 
Then this solution blows up in finite time. 

Proof. Based on the derivation in section 2 the function w(0,z,t) satisfies (35)-(36), with initial value 
w(0,z,0) = <p(z). Therefore, by Proposition 1 one concludes that for as long as the solution w(x,z,t) 
exists and is smooth must satisfies w(0,z,t) — y^-. Consequently, w(x,z,t) must lose regularity at 
sometime in the interval t € [0, 1]. 

□ 
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